Draft version July 4, 2012 

Preprint typeset using IATj^X style cmulatcapj v. 11/26/04 

gritschneder@ucolick . org 



THE ROLE OF 



OPH IN THE FORMATION AND EVOLUTION OF THE PIPE NEBULA 
FORMATION EVER ISOLATED? 



IS STAR 



(N 

o 

(N 
(N 



> 
in 
ov 

q 
o 



X 



Matthias Gritschneder 1 ' 2 and Douglas N.C. Lin 1,2 

x Kavli Institute for Astronomy and Astrophysics, Peking University, Yi He Yuan Lu 5, Hai Dian, 100871 Beijing, China; 
2 Astronomy and Astrophysics Department, University of California, Santa Cruz, CA 95064, USA 

Draft version July 4, 2012 

ABSTRACT 

We propose that the Pipe Nebula is an HII region shell swept up by the B2 IV /3 Cephei star 9 
Ophiuchi. After reviewing the morphological evidence by recent observations, we perform a series 
of analytical calculations. We use realistic HII region parameters derived with the radiative transfer 
code Cloudy from observed stellar parameters. We are able to show that the current size, mass and 
pressure of the region can be explained in this scenario. We investigate the configuration today and 
come to the conclusion that the Pipe Nebula can be best described by a three phase medium in 
pressure equilibrium. The pressure support is provided by the ionized gas and mediated by an atomic 
component to confine the cores at the observed current pressure. In the future, star formation in these 
cores is likely to be either triggered by feedback of the most massive, gravitationally bound cores as 
soon as they collapse or by the supernova explosion of 9 Ophiuchi itself. 

Subject headings: HII regions — ISM: clouds — ISM: individual objects: Pipe Nebula — ISM: structure 
— methods: analytical — methods: numerical — stars: formation 



1. INTRODUCTION 

The Pipe Nebula, first observed bv lOnishi et~aT1 (|1999f l 
is a nearby molecular clou d region. Due to its rela- 
tive pr oximity (D » 130pc, lLombardi et al.l 12006. D ~ 
145pc. lAlves fc Francoll2007t) it provides an i deal testbed 
to ob serve molecular cloud core formation (|Lada et al.1 
120081 ). The total spatial extend of the Pipe Nebular 
is roughly 14 pc x 3pc. The early measurements by 
Onis hiet al.1 (|1999f ) showed the gaseous component emit- 
ting 12 CO to have a total mass of about 10 4 M Q and the 
component emitting 13 CO to have about 3 x 10 3 M Q . In 
addition, they identified 14 core s emitting C 18 Q. More re- 
cent extinction measurements (jLombardi fc Alvesl |200H 
lLombardi et al.1 12006D increased the number of cores to 
more than 150 (lAlves et al J [20071: iMuench et al l 120071; 
iRathborne et aljl2009t iRoman-Zuhiga et al.ll2010t ). The 



total mass in the upper part of the observed area (11 pc x 
18 pc) as inferred from e xtinction measurement s is found 
to be (11000±2600)M Q (|Lombardi et al.ll2006D and thus 
in very good agreement with the previous CO observa- 
tions 1 . It is a remarkably quiescent region and there- 
fore often assumed to be the ideal case to study isolated 
star formation. Only in one tip, termed B59, obser- 
vations indicate low rate star formation and pr obably 
Jeans fragmentation (|Roman-Zuniga et al.l|2009f). The 
obser ved stars in B59 are about 2.6 Myr old (|Covev et al.l 
2010). Recent observations with the Spitzer Space Tele- 
scope confirm the ex tremely low level of s tar formation 
in the Pipe Nebula (jForbrich et all I2009D. Additional 
observation s in the n ear infrared (|R oman-Zuiiig a~et al.l 
I2009L I201dh. in NH3 ( Rathborne et all 12008(1 and X-rav 
(jForbrich et al.l l2010f) support these findings. Magnetic 
field observations (jAlvesetalJ 120081: iFrau et alj l201d: 
iFranco et alJl2010( l indicate that, especially in the Stem 



1 Denote that[Omshi et al. (1999) assume d a larger distance fo r 
the Pipe Nebula. The mass estimate of Lombardi et al. (2006) 
would be RJ 1.7 X 10 4 Mq at that distance. 



of the Pipe Nebula, the magnetic field seems to be aligned 
perpendicular to the main filament. Very recent obser- 
vations with Herschel (|Peretto et al.ll20T2T) show that in 
the Stem the sub-filaments are more grid-like, whereas 
in B59 they are more centrally condensed. 

The Pipe Nebula is located at the edge of the Sco 
OB2 association. The closest massive star is the B2 IV 
(3 Cephei star 9 Ophiuchi (HD 157056) at a projected 
distance of about 3pc from the Pipe Nebula. While 
its variability was known for a long time, only recently 
it has been target of detailed asteroseismologic studies. 
Non-rigid rotation was proven fro m multiplet observa- 
tion. Observations and models (e.g. lHandler et al"1l2005l : 
iBriquet et al.l f2007h show that 9 Oph is a triple system 
with an inner binary Stellar models (jLovekin fc Goupill 
2010) indicate that the brightest star 9 Oph A can be best 
fit by stellar models with a luminosity of log(L / L c ) = 
3.75 and an effective temperature T e g = 22590 K at an 
age igtar = 15.6 Myr. The second closest massive star is 
the B0 star r Sco at a projected distance of about 20 pc. 

Various models have been proposed to explain the for- 
mation and especially the observed core m ass function in 
the Pipe Nebula (e.g. iHeitsch et alj |2009). The effect of 
9 Oph in t riggering star forma tion has been previously 
studied by lOnishi et al.l (|1999l) . Here, they show that 
the stellar winds can not trigger star formation. How- 
ever, they did not investigate the effect of the ionizing 
radiation in the formation and the current state of the 
region. In this work, we interpret the Pipe Nebula as a 
swept up HH-region shell. In Sj2]we review the underlying 
physics. In Sj3]we present analytic models and assess the 
current state in 2) We investigate the future evolution 
in Sj5] and conclude in ££6j 

2. BASIC EQUATIONS 

In general, the evolution of an ionized region can be 
characterized by two phases. First, the ionizat ion of the 
so called Stroemgren sphere (|Stromgrenlll939t l. the vol- 
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ume which can be immediately ionized by the star before 
the system can react to the increased temperature and 
pressure. This sphere has the radius 



Rs = 



3 J, 



Ly 



1/3 



(1) 



where Jl y is the total ionizing flux of the source, no the 
number density in the surrounding and as the sum of 
the recombination coefficients for all levels besides the 
ground stage. In the second phase, the heated gas re- 
acts to its increase in temperature and therefore pres- 
sure. An approximately isothermal shock is driven into 
the surrounding medium. Under the assumption of a 
thin shock, the time evolution of the radius is given as 



(2) 



where a Si hot is the sound speed of the hot, ionized gas. 
The time evolution of the density in the ionized region is 
then given by 



(3) 



For the accumulated mass we assume that the shell con- 
tains all material swept up from the HH-region, weighted 
by the current covering fraction / of the Pipe Nebula for 
the maximum size of the bubble: 
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M(t) = — R(t)( Po - p(t))f, 



where 
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Here, ^4pi pe is the current area of the Pipe Nebular tan- 
gential to the massive star and i? m ax is maximal radius 
of the Hll-region, i.e. the current distance between the 
star and the nebula. 

A more detailed treatment has to include the fact that 
the region around 6 Oph is not spherically symmetric. 
iKrumholz fc M atzncr (2009) derived an analytic approx- 
imation for blister-type HII regions. In this case 



Rb{t) = Rs 
where ts = i?s/a S! hot, and 
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Mb(i) can then be calculated according to Eq. 2J 

3. FORMATION OF THE PIPE NEBULA 

Using these basic equations we can now derive a model 
fo r the Pipe Nebula. In the NICER extinct ion maps (Fig. 
QJ iLombardi etHI 120061 : lLada et all [2008 K a suggestive 
geometry can already be seen. The region around Oph 
(denoted by the asterisk) provides significantly less ex- 
tinction, i.e. is less dense. We interpret this as the HII 
region. 

In the following, we assume the cold gas to be at 
2coid — 10 K with a mean molecular weight of ^ co id = 
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Fig . 1. — Obser vation of th e Pipe Nebula by Lombardi et alj 
(USD, Figure 8 of lLada et all tWM) . The location of 9 Oph is 
indicated by the asterisk. The size of the circles indicates the 
pressure of the cores (in Kcm -3 ). 



1.37 ()Lombardi et al.l [ 200 6) , corresponding to a sound 
speed of a SiCO id = 0.25 km s" 1 . 2 The temperature in 
the Hll-region is resulting from the balance of heating 
by p hotoelectrons and cooling by forbidden metal lines 
(e.g. IOsterbrocklll989ft . 

We test three models A, B and C, corresponding to 
initial number densities no in the cold surrounding of 
1 x 10 3 cm~ 3 , 5 x 10 3 cm~ 3 and 1 x 10 4 cm~ 3 , respec- 
tively. In order to get a precise estimate for the radius 
and temperature of the initial Stroemgren sphere, we em- 
ploy the radiative transfer code Cloudy. Calculations 
were perfo rmed with ver s ion 8.00 of Cloudy, last de- 
scribed bv iFerland et al.l (fl99 8) . We only consider the 
most luminous star, 9 Oph A, as it is at le ast an order of 
magn itude brighter than its companions (H andler et al.l 
2005) 3 . We parametrize 6 Oph A as a black body 
with a temperature of T B ff = 22590 K and a lu minos- 
ity of log(L/L ) = 3.75 (lLovekin fc Goupill[20ll 4 . The 
code Cloudy yields R s = 0.094 pc, T hot ~ 6000 K in 
case A, R s = 0.029 pc, T hot ~ 7000 K in case B and 
R s = 0.017 pc, T hot ~ 7000 K in case C. The correspond- 
ing sound speeds are calculated assuming a molecular 
weight of phot — 0.55. 

In Fig. [2] we show the evolution in all three cases ac- 
cording to Eqn. [21 El El and [3 Solid lines correspond 
to the classical (spherical) HII region, dashed lines cor- 
respond to the blister case. The dotted horizontal line is 
the value inferred from observations. It can be directly 
seen that in all models the final radius is bigger than 
3 pc. This is especially true for th e blister case. However, 
recent numerical simulations by iGendelev Sz Krumholzl 
(2012) show that in the blister case the analytical ap- 

2 Arguably, the gas could be initially molecular, e.g. fJ. co id = 
2.35. As this just corresponds to a different choice of initial no, 
we assume it to be atomic to enable an easier comparison with the 
photo-dissociated post-shock region later on. 

3 We also neglect the influence of the B0 star r Sco, since its 
projected distance to the Pipe Nebula is 20 pc 

4 The simulations are performed under the assumption of a con- 
stant density, solar metallicity and a magnetic field of B = 10 G 
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Fig. 2. — Time evolution of the HH-region for the three different cases. Solid lines: classical (spherical) Hll-region, dashed: blister-type 
Hll-region. The lines are continued dotted once the shock has reached equilibrium with the ambient surrounding. Dotted horizontal lines: 
current day observational values. Left panel: radius of the shell. Center: swept up mass in a Pipe Nebula sized region. Right: pressure in 
the hot, ionized gas. 



proximation overestimates the radius. In their simula- 
tions it is only bigger by about 15% than the spherical 
estimate, while the analytic prediction would be 20%. 
In both cases the value is of order of the observed value 
and at the distance of the Pipe Nebula (« 130 — 146 pc) 
the discrepancy can be a projection effect, as Eqn. [5] 
and [6] give the real radius, while we can only observe the 
projected radius. 

Next, we take a look at the mass in the swept-up shell. 
We estimate the surface of the Pipe Nebula tangential to 
the direction of the star to be the same as the projected 
size of the Nebula, i.e. ^4pi pe = 14 pc x 3pc. This is of 
course ad hoc and only an order of magnitude estimate. 
However, panel 2 of Fig. [2] shows clearly that in the HII- 
region scenario it is no problem to accumulate enough 
mass to be consistent with the observed value in the Pipe 
Nebula (M » 2.2 x 10 3 M , see ©. 

Another constraint for the models is the pressure re - 
quircd to confine the observed cores. lLada et al.l (|2008l ) 
estimate the required average pressure at the surface of 
the cores to be ~ 8 x 10 4 Kcm -3 . In our model, this 
pressure has to be initially supplied by the ionized gas 
(for a detailed discussion see . From Fig. [3] it can be 
seen that it is challenging to maintain this pressure in 
our model for the entire time-span of 15.6 Myr. Never- 
theless, the analytical value is close to the observed value 
in all three cases. In addition, a closer investigation of 
Fig. Q] shows, that the cores in the swept up shell have 
lower pressures, whereas the cores in the central region 
have higher pressure. We discuss this further in 21 

The shock front will stagnate as soon as the pressure 
in the neutral gas equals the pressure in the ionized gas. 
This is reached in Cases B and C. Once the front stag- 
nates, the corresponding lines in Fig. [2] are continued 
dotted. All three quantities are remarkably close to the 
observations at this point. This leads to a straightfor- 
ward interpretation. The HII region expanded until stag- 
nation. At the same time, part of the shocked gas and 
the post-shock gas expand, as they are preheated by the 
photodissociating photons of the B-star to about 100 K 
(see ID). The increased motion (a SjWarm = 0.78kms _1 ) 
leads to a pile-up of material towards the solid boundary 
of the ionized gas and to a broadening of the shock in the 
opposite direction. Thus, the shock broadens with half 



is, warm, reaching about 3.2 pc in the 8 Myr since stag- 
nation in Case B. The density at stagnation in the cold 
shock is ng and will drop according to the increase in tem- 
perature. In case B, with a new temperature of 100 K, 
this corresponds to 500cm -3 . Taking into account the 
additional pile up this is very close to the value inferred 
from the observations (^atomic = 774cm -3 , see below). 

In addition, the equations in <J3] are derived under the 
approximation of a thin shock. This approximation will 
fail eventually. Then, the speed of the ionization front 
will be lower than the analytical estimate. Therefore, 
the radius and swept up mass are best viewed as an up- 
per limit. Correspondingly, the pressure is a lower limit, 
since a smaller region will lead to higher densities and 
thus to a higher pressure in the ionized gas. 

4. CURRENT STATE 

We now continue to discuss the current state of the 
Pipe Nebula in more detail. The most obvious con- 
straints are the mass, pressure and current size. 

To asses the current mass, it is necessary to take a 
detailed look at the observations. The total mass in 
the upper part of the observed region (b > +3°) is 
inf erred to be 10 4 M[T) fr om extinction measurements 
by lLombardi et al.l (|2006t ). They provide a detailed 
comparison with t he previous 12 CO-measurements by 
lOnishi et al.l ()1999t ) and find encouraging agreement. To 
estimate the mass inside the Pipe Nebula itself it is first 
necessary to look at the ma s s in t he dense component 
alone. From lLombardi et al.l (|2006h Fig. 7, it is reason- 
able to assume that the dense component has an extinc- 
tion Ak > 0.4. Inside this isoex tinction contour there 
is about 40% of the total mass (jLombardi et al.l 120061 
their Fig. 27). About half of this mass is in the Pipe 
Nebula, therefore the denser component of the nebula 
itself has a mass of about 2.2 x 10 3 M Q . This is in agree- 
ment with th e 3x1 3 Mq in the entire region inferred by 
lOnishi et al.l (|1999f ) from 13 CO, which traces the denser 
gas (e.g. their Fig. 5). The to tal mass in the cores is 
about 250 M (|Lada et al.ll2008l their Table 1). 

From these masses, we can derive the structure of 
the Nebula. The cores have a known pressure of w 
8 x 10 4 Kcm -3 . The ionized gas can supply this pres- 
sure. Independent of the uncertainties in <J3J we per- 
form a Cloudy simulation on the current static state. It 
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shows, that 9 Oph A can ionize gas at no = 8 cm 3 up 
to a radius of 2.75 pc to T = 5000 K. With u ion = 0.55 
the pressure is then 7.3 x 10 4 Kcm~ 3 , which is a bit less 
than the observed average density of the cores. However, 
a closer look at Fig. [1] reveals, that the cores along the 
rim generally have a lower pressure (i.e. smaller circles), 
whereas the high pressure cores are towards the central 
region. A straightforward assumption is that the atomic 
component in between the cores and the ionized gas is in 
pressure equilibrium with both, acting as a buffer. As- 
suming a volu me of 14 pc x 3 pc x 2 pc for the Pipe Nebula, 
and (J, = 1.37 (jLombardi et al.ll2006l ). this corresponds to 
^atomic = 774 cm~ 3 . To be in pressure equilibrium, the 
atomic phase has to be at T = 100K. Thus, the bulk 
of the mass (« 87.5%) is in the atomic phase. To keep 
the medium at this temperature, significant heating is 
required. The heating can be suppli e d by the diffuse ra- 
diation. In addition, lArthur et al.l f|2011h showed that 
for B-stars the photo-dissociation-region (PDR) is wider 
than for O-stars. This is due to the peak at lower ener- 
gies in the spectrum. Therefore, in B-stars, the ratio of 
photo-dissociating photons to ionizing photons is higher. 
As a consequence, the photo-dissociation, i.e. the disrup- 
tion of H2 molecules, outruns the ionization front, leading 
to a warmer region at about 100 K with a thickness of 
about 30% of the radius of the HII region. 

Concerning the high pressure cores in the central re- 
gion, it is probably best to assume that this region is 
currently undergoing gravitational collapse. Possibly, the 
swept up HH-shell encountered a denser region and the 
resulting density enhancement lead to a gravitational un- 
stable region. Thus, the cores in the region are gravita- 
tionally bound. 

The width of the Pipe Nebula can be explained by 
assuming the shock and post-shock region expand after 
stagnation with ^a SiWar m (see $3]). Therefore, the Case C 
would lead to 4.4pc and Case B to 3.2 pc, as observed. 

5. FURTHER EVOLUTION AND IMPLICATIONS FOR THE 

IMF 

As the bulk of the mass is in the atomic phase, this 
phase will most likely determine the future evolution of 
the Pipe Nebula. A possible scenario is triggered star 
formation from within. When the most massive clumps 
turn into stars, they will ionize and thereby heat the 
atomic phase. This will in turn increase the pressure on 
the cold cores, forcing them to form stars. 

In the context of the transition from the Core Mass 
Function (CMF) to the Initial Mass Function (IMF) it 
is very interesting to look at the increase of the ambient 
pressure quantitatively. If the radiation of the newborn 
stars increases the background pressure by a factor of 
two, the shift from the I MF to the CMF could be ex- 
plained (|Lada et al.ll2008t) . 

There are three cores of about 10 M© and one with 
« 20 M Q in the Pipe Nebula. Assuming that the entire 
most massive cores end up in a single stars, the further 
evolution will be similar to our model A. As can be di- 
rectly seen from Fig. El there would be a shock pro- 
ceeding through the Pipe nebula. However, due to the 
number density in the atomic gas, the evolution of this 
shock would take several Myr, whereas star formation in 
general seems to be more synchronized. For the 10 M© 
cores, the evolution would be even slower. Furthermore, 



this triggering scenario would leave the top end of the 
IMF unaffected as these cores turn into stars before they 
themselves increase the ambient pressure. 

A more violent scenario will occur if 9 Oph explodes 
in a supernova before the massive cores turn into stars. 
Given the relatively small distance, the Pipe Nebula will 
be hit by a blast wave in the Sedov phase within a 
few kyr after the explosion. Mostly likely, the smaller 
cores will get disrupted, while the bigger cores will get 
triggered into collapse w ithin another few kyr (see e.g. 
iGritschneder et aLll2012[) . This could indeed explain the 
small age spread observed in star forming regions. 

Another environmental effect is the influence of the 
Sco OB2 asso ciation. As pointed out pr eviously (e.g. 
Ide Geusl 119921 ) and discussed recently by iPeretto et al.l 
(2012), the star forming region termed B59 is pointing 
towards that association and resembles a bow shock like 
structure. Howeve r, the age estimate for the young stars 
in B59 of 2.6 Myr (jCovev et all 120101 ) suggests that this 
interaction triggered the stars in the past and is not cur- 
rently driving star formation further down the stem. 

6. DISCUSSION & CONCLUSIONS 

There are several uncertainties involved in the model. 
The main uncertainty lies in the assumption of a thin 
shock for the analytic solution. Once the shock thickens, 
it will lead to a smaller region and a less dense shell with a 
higher pressure, as discussed above. Another possibility 
to constrain the H II region to the curre nt size would 
be magnetic fields. iKrumholz et al.l ()2007l ) were able to 
demonstrate that the evolution of a shell stagnates earlier 
for a magnetized medium since the Alfven speed is higher 
than the sound speed (see also lArthur et al]l2011| ). 

Another uncertainty lies in the assumed initial condi- 
tions for the Cloudy models, e.g. the assumed abun- 
dances, grains and magnetic fields will influence the out- 
come. However, for the temperature and the size of the 
HII region, these simulations give more realistic answers 
than the simple assumptions leading to Eq. [2] Further- 
more, the most important input parameter - besides the 
stellar parameters taken directly from observations - is 
still the density, as it enters with Uq. 

The assumed time scale for the sweeping up of the 
Pipe Nebula is « 15 Myr. This is comparable to the free- 
fall time of a cloud so a remaining question is why the 
nebula did no t collapse. One poss ible solution might be 
the finding of I Arthur et al.l (|2011| ) that B-stars have an 
extended PDR. In this broad region behind the ionization 
front the temperature is c lose to 100K. As e.g. shown by 
IGritschneder et all (|2010l ). a higher temperature in the 
cold gas is prohibitive of structure formation and thus 
hindering collapse. 

Besides these uncertainties, the scenario of 9 Oph swip- 
ing up the Pipe Nebula presented in this paper can, es- 
pecially in Case B, successfully explain: 

1 . The observed morphology of the Pipe Nebula, espe- 
cially the spheroidal shape of the less dense region. 

2. The current width, as a shock broadening with 
5Qs S)Warni = 0.44 km s -1 will reach a width of 3.2 pc 
in the 8 Myr since the stagnation. 

3. The current mass and size can be easily reached. 
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4. More importantly, the pressure to confine the cores 
can be supplied. 

Especially the pressure of the cores is otherwise puzzling. 
Up to now, the only possible explanation for this confine- 
ment was the self-gravity of the cloud. This is highly un- 
likely, as the cores in a self-gravitating system are the first 
instances to react to the collapse and therefore should be 
bound, whereas most of the cores are observed to be un- 
bound. 



7. ACKNOWLEDGEMENTS 

M.G. acknowledges funding by the Alexander von 
Humboldt Foundation in form of a Feodor-Lynen Fel- 
lowship and by the China National Postdoc Fund Grant 
No. 20100470108 and the National Science Foundation 
of China Grant No. 11003001. D.N.C.L. acknowledges 
funding by the NASA grant NNX08AL41G. 



REFERENCES 



Alves, F. O., & Franco, G. A. P. 2007, A&A, 470, 597 
Alvcs, F. O., Franco, G. A. P., & Girart, J. M. 2008, A&A, 486, 
L13 

Alves, J., Lombardi, M., & Lada, C. J. 2007, A&A, 462, L17 
Arthur, S. J., Henney, W. J., Mellema, G., de Colle, F., & Vazquez- 

Semadeni, E. 2011, MNRAS, 414, 1747 
Briquet, M., Morel, T., Thoul, A., Scuflaire, R., Miglio, A., 

Montalban, J., Dupret, M.-A., & Acrts, C. 2007, MNRAS, 381, 

1482 

Covey, K. R., Lada, C. J., Roman-Zufuga, C, Muench, A. A., 

Forbrich, J., & Ascenso, J. 2010, ApJ, 722, 971 
de Geus, E. J. 1992, A&A, 262, 258 

Ferland, G. J., Korista, K. T., Verner, D. A., Ferguson, J. W., 
Kingdon, J. B., & Verner, E. M. 1998, PASP, 110, 761 

Forbrich, J., Lada, C. J., Muench, A. A., Alves, J., & Lombardi, 
M. 2009, ApJ, 704, 292 

Forbrich, J., Posselt, B., Covey, K. R., & Lada, C. J. 2010, ApJ, 
719, 691 

Franco, G. A. P., Alves, F. O., & Girart, J. M. 2010, ApJ, 723, 146 
Frau, P. et al. 2010, ApJ, 723, 1665 
Gendelev, L., & Krumholz, M. R. 2012, ApJ, 745, 158 
Gritschneder, M., Burkcrt, A., Naab, T., & Walch, S. 2010, ApJ, 
723, 971 

Gritschneder, M., Lin, D. N. C, Murray, S. D., Yin, Q.-Z., & Gong, 

M.-N. 2012, ApJ, 745, 22 
Handler, G., Shobbrook, R. R., & Mokgwetsi, T. 2005, MNRAS, 

362, 612 



Hcitsch, F., Ballesteros-Paredes, J., & Hartmann, L. 2009, ApJ, 
704, 1735 

Krumholz, M. R., & Matzner, C. D. 2009, ApJ, 703, 1352 
Krumholz, M. R., Stone, J. M., & Gardiner, T. A. 2007, ApJ, 671, 
518 

Lada, C. J., Muench, A. A., Rathbornc, J., Alves, J. F., & 

Lombardi, M. 2008, ApJ, 672, 410 
Lombardi, M., & Alves, J. 2001, A&A, 377, 1023 
Lombardi, M., Alves, J., & Lada, C. J. 2006, A&A, 454, 781 
Lovekin, C. C, & Goupil, M.-J. 2010, A&A, 515, A58 
Muench, A. A., Lada, C. J., Rathborne, J. M., Alves, J. F., & 

Lombardi, M. 2007, ApJ, 671, 1820 
Onishi, T. et al. 1999, PASJ, 51, 871 

Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae and active 

galactic nuclei (University Science Books) 
Peretto, N. et al. 2012, A&A, 541, A63 

Rathborne, J. M., Lada, C. J., Muench, A. A., Alves, J. F., 

Kainulainen, J., & Lombardi, M. 2009, ApJ, 699, 742 
Rathborne, J. M., Lada, C. J., Muench, A. A., Alvcs, J. F., & 

Lombardi, M. 2008, ApJS, 174, 396 
Roman-Zuniga, C. G., Alves, J. F., Lada, C. J., & Lombardi, M. 

2010, ApJ, 725, 2232 
Roman-Zuniga, C. G., Lada, C. J., & Alves, J. F. 2009, ApJ, 704, 

183 

Stromgren, B. 1939, ApJ, 89, 526 



